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The KamLAND and Borexino experiments have observed, each at ~ 4(t level, signals of electron 
antineutrinos produced in the decay chains of thorium and uranium in the Earth's crust and mantle 
(Th and U geoneutrinos). Various pieces of geochemical and geophysical information allow an 
estimation of the crustal geoneutrino flux components with relatively small uncertainties. The 
mantle component may then be inferred by subtracting the estimated crustal flux from the measured 
total flux. To this purpose, we analyze in detail the experimental Th and U geoneutrino event rates 
in KamLAND and Borexino, including neutrino oscillation effects. We estimate the crustal flux at 
the two detector sites, using state-of-the-art information about the Th and U distribution on global 
and local scales. We find that crust-subtracted signals show hints of a residual mantle component, 
emerging at ~ 2.4(t level by combining the KamLAND and Borexino data. The inferred mantle 
flux slightly favors scenarios with relatively high Th and U abundances, within ±fcr uncertainties 
comparable to the spread of predictions from recent mantle models. 



I. INTRODUCTION 

The decay chains of uranium (U), thorium (Th), and potassium (K) in the Earth's interior provide intense sources 
of terrestrial heat and, at the same time, of low-energy electron antineutrinos (i7g) — the so-called geoneutrinos [ij. 
Geoneutrinos from Th and U (but not from K) decay are detectable via the inverse beta decay (IBD) reaction, 

i7e-|-p^n-|-e+ (£;j. > 1.806 MeV) , (1) 

and have recently been observed at ~ 4cr level both in the KamLAND (KL) [31 and in the Borexino (BX) Q 
experiments. The KL and BX measurements represent first steps of a long-term research program vifhich, by bridging 
particle physics and Earth science, vifill provide unique clues on fundamental geophysical and geochemical issues jJ-Q- 

Indeed, the geoneutrino flux and its spectrum encode relevant information about the distribution of radiogenic 
elements in the crust and in the mantle, which are thought to be the main Th and U reservoirs ^} In particular, 
the total Fg flux probes the total amount of radiogenic elements in the Earth, while the energy spectrum is sensitive 
to the different Th and U components [ij. In principle, the angular spectrum (not yet experimentally accessible fzl) 
can probe the different mantle and crust source geometry [8]. Repeating such measurements at different locations 
can thus help to distinguish the site-dependent crustal components from the (approximately) site-independent mantle 
component of the flux, which can also be more directly probed at oceanic sites [9|]. 

Extracting such information is not straightforward, since the geoneutrino flux represents a volume integral over 
Th and U abundances, weighted by the inverse square distance, and modulated by the IBD cross section and j7g 
oscillation probability (see [1] for details). While the latter two ingredients are known with good accuracy, the volume 
distribution of Th and U is subject to relatively large uncertainties, especially in the mantle Therefore, in order to 
disentangle interesting pieces of information from particle physics data (geoneutrinos) , one needs an interdisciplinary 
approach, including supplementary constraints or assumptions from Earth science (geophysics and geochemistry). 
For instance, in order to constrain the radiogenic heat flux, one may exploit its (partly mo del- dependent) covariances 
with the total geoneutrino flux [H, oi, in order to probe meteoritic expectations for the Th/U abundance ratio 
[ll| . one may assume that KL and BX experiments probe the same average Th/U in first approximation [12]. 

In this work we apply such an interdisciplinary approach to infer the mantle component of the geoneutrino flux, 
which we obtain by subtracting accurately estimated crust components from the total measured fluxes. In particular, 
concerning particle physics data, we perform a detailed analysis of the total Th and U geoneutrino fluxes measured 
in KL and BX, including oscillation effects (Sec. II). Concerning Earth science data, we estimate the different crustal 
flux components in the two experiments, using state-of-the-art geochemical and geophysical information about the 



Geochemical arguments disfavor significant amounts of Th and U in the Earth's core, see [3| and refs. therein. If the core were 
(hypothetically) a geoneutrino source, the "mantle" fluxes estimated in this work should be interepreted as "core+mantle" fluxes. 
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crust, on both global and local scales (Sec. III). The mantle component in KL and Borexino is then obtained by 
subtraction (mantle — total — crust) . Within the reasonable assumption of site- independent mantle flux, the KL and 
BX results can be combined, yielding a mantle signal at > 2a (Sec. IV). In comparison with current models of the 
mantle, the signal best fit favors scenarios with relatively high Th and U mantle abundances, within ±lcr uncertainties 
comparable with the spread of model predictions (Sec. V). The main results are summarized in Sec. VI. Statistical 
details and side results of our analysis are confined to Appendix A and B, respectively. 

II. PARTICLE PHYSICS INPUT: ANALYSIS OF KL AND BX DATA 

In this Section we update the KL and BX geoneutrino data analysis discussed in [l^ , by including recent constraints 
from world neutrino oscillation data as in [l^ , and the latest KL geoneutrino data from . With respect to , the 
BX data 0] are the same, while the updated KL data are significantly more accurate. We use the KL energy spectra 
and detection efficiencies for 77e events as shown in Fig. I of [2] , within the same statistical approach discussed in [l^l ■ 
In particular, the fit to KL and BX data involves a 7-dimensional manifold. 

Parameters = {W, 0i2, ^is; i?(Th)KL, R{V)kl, i?(Th)Bx, R{V)bx} ■ (2) 

where the four i?'s represent the KL and BX event rates from Th and U geoneutrinos, expressed in Terrestrial 
Neutrino Units (1 TNU — I0~^^ events per target proton per year). It is useful to remind that a reference (oscillated) 
Ve flux (j) = lO^/cmVs generates i?(Th) = 4.04 TNU and i?(U) = 12.8 TNU, while a mass abundance ratio Th/U = 
a(Th)/a(U) corresponds to i?(Th)/i?(U) = 0.0696 • Th/U for a homogeneous source [3. 

In Eq. ([2]), the mass-mixing oscillation parameters (Sm^, 9i2, ^13) govern the flavor survival probability Pee of both 
geo-Fe and background reactor V^. , 

Pee = P{Ve ^ Ve) = cos^ (^1 - sin^ 2^12 sin^ (^^^)) + siu* ^13 , (3) 

L and E being the Ve path length and energy, respectively, in natural units. As in we make the reasonable 
approximation of oscillation-averaged Pee for geoneutrinos, 

{Pee) ^ COS^ 013 (^1 " ^ siu^ 2^12^ + siu^ ^13 . (4) 

We have verified a posteriori that, within current uncertainties, removing this simplifying assumptions does not spoil 
our main results; see Appendix A for details. 

For (Pee), we adopt the reference la ranges sin^ 612 ~ 0.306 ± 0.017 and sin^ 6*13 ~ 0.021 ± 0.007 from the global 
analysis of oscillation data (from solar, atmospheric, accelerator, and reactor v experiments) in [l3| . implying: 

(Pee) =i 0.551 ± 0.015 (Icr) . (5) 

The 3% fractional uncertainty of (Pee) is much smaller than other (mainly experimental) errors affecting current 
geoneutrino event rates. In any case, it is properly taken into account by marginalization in the global At of 
oscillation plus geoneutrino data, which yields measured event rates R which consistently include all the uncertainties 
from particle physics inputs.^ 

Figure 1 shows the results of our analysis of the total (Th-|-U) geoneutrino event rate R in KL and BX, in terms 
of standard deviations (N^ — \J Ax^) from the best fit. The null hypothesis of vanishing geoneutrino flux (i? = 0) is 
rejected at ~ 4.2 by both KL and BX, in good agreement with the corresponding official results in 0] and Q; see 
also Appendix B. Note that the KL and BX curves are not linear and symmetric, as it would be the case for gaussian 
errors: indeed, the asymmetry and nonlinearity increase for decreasing event rate R (especially for the low-statistics 
BX data), as a result of Poisson fluctuations, whose statistics is properly accounted for in our analysis [l2|. As a 
consequence, the current KL and BX event rates cannot be simply summarized in terms of central values and ztltr 
errors: non-gaussian uncertainties must be properly taken into account, in order to exploit the full potential of the 
available geoneutrino data. 



^ After our work was completed, high-precision measurements of sin^ 613 were announced by two reactor neutrino experiments: Daya Bay 
(sin^ 613 = 0.024 ±0.004) ,141] and RENO (sin^ 6»i3 = 0.029 ±0.006) [ll. The currently recommended \a range from world reactor data, 
sin^ 013 ~ 0.025 ± 0.003 jl^, is fully contained in the Icr range used herein: sin^ 013 ~ 0.021 ± 0.007 [l3|. An updated global analysis of 
the neutrino mass-mixing parameters, including Daya Bay, RENO and other recent oscillation data, is reported elsewhere I?]. In any 
case, such update would have negligible effects on the geoneutrino results presented in this work. 
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R(Th + U) [TNU] 

FIG. 1: Constraints on the total geoneutrino event rate in KL and BX, in terms of standard deviations from the best fit. 



In order not to loose precious information, in the following we shall mostly refer to the full data analysis in terms 
of separate (not summed) Th and U event rates. Figure 2 shows, in particular, the main results of our analysis of 
KL and BX data, in the plane charted by the geoneutrino event rates -R(U) and i?(Th) for KL (upper panel) and BX 
(lower panel). In both panels, the thick dot represents the best fit, while the curves correspond to the na contours 
(Ax^ = n^, for n = 1, 2 and 3), whose projections provide the na bounds for the corresponding parameter [l^. Note 
that, in Fig. 2, the la contour for KL is closed, i.e., KL can separate the Th and U components at the (weak) level of 
~ 1(7, due to higher statistics and better spectral information. However, already at 2ct, the KL and the BX contours 
are no longer closed; indeed, the strong anticorrelation of the na isolines reflects the fact that the KL and BX spectra 
are currently more sensitive to the total Th+U flux than to their separate Th and U components which, to some 
extent, can be traded one for the other. Figure 2 also anticipates schematically the subtraction of the estimated 
crustal rates (with their associated ±3cr errors), as discussed in detail in the next Sections. 
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FIG. 2: Results of our analysis of geoneutrino rates _R(Th) and -R(U) expressed in TNU, for both KL (upper panel) and BX 
(lower panel). The curves represent the na contours (Ax^ — n^) around the best fit (thick dot). Also shown is the shift of the 
best fit points after subtraction of the estimated crustal components (with their ±3(t errors). See the text for details. 



4 



III. EARTH SCIENCE INPUT: CRUSTAL FLUX ESTIMATES IN KL AND BX 

In order to estimate the crustal geoneutrino flux we need a global model for the Earth crust and a sufficiently 
detailed model for the local contribution. Indeed, due to the inverse square law for the flux, the crust portions within 
and outside a radius of 0(500) km from the detector provide comparable flux contributions in both KL and BX 
In particular, for both KL and BX, we use an accurate description of the local crust extending over ^2.5 x 10^ km^ 
and down to 30-40 km of depth (Moho surface), which contributes to ^40% of the total geo-v signal. For farther 
portions of the crust, a coarser description in terms of 2° x 2° tiles is sufficient. 

In this section we discuss global and local properties of the crust, building upon previous works on the subject 
[H, particular, we report our calculation of the crustal Th and U fluxes at KL and BX, with estimated la 

uncertainties of 0(10%). It should be noted, however, that precise estimates for such uncertainties play no significant 
role in this work: even if the crustal flux errors were all conservatively doubled or tripled, the final results for the 
mantle signal would only change by a tiny fraction of one standard deviation (see Appendix A) . 

A. Global model of the crust 

Our global model for the crust is based on a geophysical 2° x 2° tiled map [1^, as well as on the Th and U mass 
abundances recommended in (2l| for sedimentary layers and in [2^ for the upper, middle, and lower crust. For 
the lower crust, the values in the literature encompass a large interval: we adopt a mean value together with an 
uncertainty indicative of the spread of published values [ij . After removal of a tiny portion of crust extending for a 
few hundred km around each site (local or "LOG" contribution, as defined in the next subsections), the global model 
is used to evaluate the remaining geoneutrino flux at both KL and BX (rest-of-the-crust or "ROC" contribution). 

Table I reports the input total mass and adopted values and ±lcr errors for radiogenic abundances in each global 
reservoir, as well as the output ROC event rates in KL and BX, multiplied by the central value of the probability in 
Eq. (O. [The Pee errors are already accounted for in the experimental data fits of Figs. 1 and 2.] Total errors in the 
last row are obtained by summing partial errors in quadrature. 

TABLE I: Inputs and outputs of the global model of the crust adopted in this work. The first four columns report, for each 
reservoir, its mass M and the adopted Th and U mass abundances. The last four columns report the estimated rest-of-the-crust 
(ROC) event rates, as obtained by excluding from the total crust the local (LOG) portions defined in the text. Quoted errors 
are at Icr. 



Reservoir 


M [10^2 kg] 


a(Th) [/ig/g] 


«(U) [/ig/g] 


ROC rates for KL 
i?(Th) [TNU] i?(U) [TNU] 


ROC rates for BX 
i?(Th) [TNU] i?(U) [TNU] 


Sediments 


0.11 


6.9 ±0.8 


1.7 ±0.2 


0.10 ±0.01 


0.34 ± 0.04 


0.22 ±0.03 


0.82 ±0.09 


Upper crust 


0.70 


10.5 ± 1.0 


2.7 ±0.6 


0.99 ±0.10 


3.64 ± 0.80 


1.66 ±0.16 


6.42 ± 1.43 


Middle crust 


0.71 


6.5 ±0.5 


1.3 ±0.4 


0.62 ± 0.05 


1.80 ±0.56 


1.11 ±0.09 


3.32 ± 1.02 


Lower crust 


0.66 


3.7 ±2.4 


0.6 ±0.4 


0.34 ±0.22 


0.80 ±0.54 


0.59 ±0.39 


1.44 ± 0.96 


Oceanic crust" 


0.60 


0.22 ±0.07 


0.10 ±0.03 


0.02 ±0.01 


0.11 ±0.04 


0.01 ± .003 


0.07 ± 0.02 


Total 








2.07 ±0.25 


6.71 ± 1.12 


3.72 ± 0.43 


12.07 ±2.00 



"For the oceanic crust, uncertainties are taken from private communication with R. Rudnick. 



From the global model we estimate the following contributions to the radiogenic crustal heat: i/(Th) = 4.02 ± 
0.47 TW and H{\J) = 3.40±0.56 TW. Assuming a mass ratio K/U - 13000 [H], the additional K contribution would 
be iJ(K) ~ 1.5±0.2 TW and the total estimated crustal heat would amount to i7(Th + U±K) ~ 8.9±1.2 TW, where 
errors have been added linearly, due to the high positive correlations among the three radiogenic element abundances. 

B. Local model of the crust around KL 

The local crust at the KL site (Kamioka) is defined in terms of six 2° x 2° tiles, supplemented with geochemical 
information on a 0.25° x 0.25° grid and on a detailed map of the crust depth [l^ . The possible effects of the subducting 
slab beneath Japan are considered, and the uncertainties arising from the debated (continental or oceanic) nature of 
the crust below the Japan Sea are also taken into account. The maximal and minimal excursions of various inputs 
and uncertainties are taken as a proxy for the ±3(7 error range. 

Table II summarizes the estimated LOC contributions to the geoneutrino signal in KL, together with their estimated 
la errors; total errors are obtained by summing in quadrature. Further details can be found in [ll. [isj. 
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TABLE II: Local (LOG) contributions to the geoneutrino signal in KL. Quoted errors are at la. 



Reservoir 


i?(Tli) [TNU] 


i?(U) [TNU] 


Six tiles" 


3.20 ±0.37 


11.17 ±0.65 


Subducting slab 


0.90 ±0.27 


2.02 ±0.61 


Japan sea 


0.09 ± 0.03 


0.34 ±0.10 


LOG total 


4.19 ±0.46 


13.53 ±0.90 



"The six-tiles errors include uncertainties on the crust composition, depth, and map discretization. 



C. Local model of the crust around BX 

The local crust at the BX site (Gran Sasso) is defined in terrns of a 2° x 2° central tile (CT) and of the rest of the 
region (RR) formed by the surrounding six tiles minus the CT jl9| . Geophysical features of the local crust (geometry, 
density, seismic velocities etc.) are reported in [l^. From a geochemical viewpoint, the CT sedimentary is a mixture 
of four main reservoirs, which have been probed by direct sampling of Th and U abundances. In the upper and 
lower crust one can recognize two components (felsic and mafic rocks) which are also probed by direct measurements. 
Average elemental abundances for the two groups were calculated and seismic arguments used in order to fix their 
relative amounts within the upper crust. In the lower crust, the fraction of felsic and mafic rocks was estimated on 
the basis of geophysical and geochemical information. The same Th and U abundances were assumed in the CT and 
in the RR. The maximal and minimal excursions of various input values and uncertainties were taken as a proxy for 
the ±3(T error range. 

Table III summarizes the estimated LOG contributions to the geo-v signal in BX, together with their la errors; 
total errors are obtained by summing in quadrature. Adopted abundances are also reported. See [l^ for details. 



TABLE III: Local (LOG) abundances and contributions to the geoneutrino signal in BX. Quoted errors are at la. 



Reservoir 


a(Th) [fig/g] 


a(U) [Mg/g] 


R{Th) [TNU] 


J?(U) [TNU] 


Sediments 


2.00 ±0.17 


0.80 ±0.07 


0.40 ± 0.04 


2.53 ±0.21 


Upper crust 


8.1 ± 1.6 


2.20 ±0.43 


1.21 ±0.24 


4.94 ± 0.96 


Lower crust 


2.6 ± 1.2 


0.30 ±0.10 


0.25 ±0.11 


0.34 ±0.11 


LOG total 






1.86 ±0.27 


7.81 ±0.99 



D. Estimated crustal geoneutrino rates and uncertainties at KL and BX 

Table IV summarizes the results of this section, in terms of LOG, ROC and total event rates in KL and BX, 
with errors added in quadrature. The fractional uncertainties, all of 0(10%), are not crucial in the context of our 
analysis, which is dominated by experimental errors (see below and Appendix A). In the future, however, it might 
be useful to address some geophysical uncertainties — such as those related to the global crust thickness — which have 
been neglected herein, in comparison with the larger geochemical abundance uncertainties. 

Concerning the central values, it should be noted that the total estimated rates in BX are lower than in KL, 
contrary to previous estimates and expectations (based on the fact that BX is surrounded by thicker crust than KL) 
l|, p] . This counterintuitive result is mainly driven by the improved description of the local BX crust performed in 
19|, leading to a significantly deeper sediment layer, and to a depletion of Th and U in all local crust reservoirs, as 
compared to previous estimates using no (or less accurate) local models. 



TABLE IV: Summary of LOG, ROG and total crust contributions to the geo-i/ signal in KL and BX. Quoted errors are at la. 





KL event rates 




BX event rates 




Reservoir 


i?(Th) [TNU] 


R(U) [TNU] 


i?(Th) [TNU] 


i?(U) [TNU] 


LOG 


4.19 ± 0.46 


13.53 ±0.90 


1.86 ±0.27 


7.81 ±0.99 


ROG 


2.07 ±0.25 


6.71 ± 1.12 


3.60 ± 0.43 


12.07 ±2.00 


Grust total 


6.26 ± 0.52 


20.24 ± 1.43 


5.46 ±0.51 


19.88 ± 2.24 
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FIG. 3: Mantle geoneutrino rates obtained from crustal rate subtraction in KL (upper panel) and BX (middle panel), as well as 
in their combination KL+BX (lower panel). The curves represent the na contours around the best fit. In the lower panel, lines 
of constant Th/U abundance ratio are also shown. The null hypothesis (no mantle signal) is disfavored at > 2a by KL+BX. 



IV. CRUST SUBTRACTION AND GEONEUTRINO RATES FROM THE MANTLE 

In the previous two sections we have determined, for both KL and BX, the experimental total rates R (Fig. 1) 
and the theoretical crustal rates i?crust (Table IV) due to Th and U geoneutrinos, together with the associated 
uncertainties. In this section we infer the mantle rates, by subtracting the estimated crustal components from the 
experimental total rates, 

_R(mantle) = _R(total, exp) — i?(crust, theo) . (6) 

Going back to Fig. 2, in each panel, crustal rate subtraction is graphically shown as a shift of the best-fit point by a 
vector (ending with crossed error bars) defined by the total Th and U crust rates and their ±3ct errors (see Table IV). 
In principle, such estimated crust uncertainties should be properly combined with the experimental rate uncertainties. 
However, the latter are currently larger by (more than) an order of magnitude, and dominate the error associated to 
the subtraction procedure. Indeed, in Appendix A we demonstrate that, in the present context, crustal rate errors 
are practically insignificant, even if they are conservatively inflated by a factor of a few. As a consequence, the na 
contours for the mantle rates can be simply obtained by rigidly shifting the na contours in Fig. 1, together with the 
best fit point, along the slanted line. Numerically, this amounts to a simple translation of the two-dimensional 
functions for both KL and BX. 

Figure 3 shows the allowed na regions of the mantle rates resulting from crust subtraction in KL (upper panel) and 
in BX (middle panel). In both cases, positive values for the Th and U mantle rates are preferred. The fact that the 
best fit rates are well within the physical region (and none of them is negative) is an encouraging sanity check of the 
crust subtraction procedure. More precisely, we find that the null hypothesis of no mantle signal (i.e., the origin of 
the axes) is disfavored at about 1.7a level in KL and 2.0a in BX. Moreover, the allowed regions in KL and BX largely 
overlap. Therefore, we have obtained two consistent hints at > 1.7a level in favor of a geoneutrino signal coming from 
the mantle. Nonzero mantle fluxes were also suggested (but not quantified) by the analyses in 0, ■ 



7 



The above hints can now be properly combined under the assumption of site-independent mantle flux, which is 
justified in the absence of significant indications in favor of local mantle anomalies below the KL or BX sites. (Possible 
geochemical mantle anomalies are still debated and, in any case, they are strongly model dependent [ISl-) Under 
such hypothesis, the combination of KL and BX constraints on mantle geoneutrino rates amount to summing the 
corresponding functions. The results are shown in the lower panel of Fig. 3. The mantle signal common to the 
KL and BX emerges now with greater statistical significance, the null hypothesis being rejected at 2.4cr (about 98.4% 
C.L.). This result represents an encouraging first step towards a better understanding of the mantle via geoneutrinos, 
and can already provide valuable indications in comparison with various mantle models, as shown in the next Section. 

Of course, the inferred mantle signal is not yet accurate enough to probe more detailed issues, such as the mantle 
Th/U ratio. In particular, in Fig. 3, the KL+BX contour at la appears to be compatible with any possible Th/U 
mantle ratio (some isolines being shown to guide the eye), so that the preference for Th/U ~ 8 is not statistically sig- 
nificant. A future reduction of the experimental errors (which is conceivable with longer exposures, better background 
rejection, and additional experiments) would be desirable to get a mantle geoneutrino signal with greater impact for 
geophysics and geochemistry. 



COMPARISON WITH MANTLE MODELS 



In this Section we compare the inferred mantle geoneutrino signal with predictions derived from various published 
mantle models, hereafter referred to as: Turcotte and Schubert 2002 \25], Anderson 2007 Palme & O'Neill 2003 



27 1, Allegre et al. 1995 [28|, McDonough & Sun 1995 [29|, Lyubetskaya & Korenaga 2007 [3d], Javoy et al. 2010 
3l|. It should be remarked that, in general, the models mainly focus on the "primitive mantle" (PM, with mass 
MpM — 4.03 X 10^'* kg) prior to differentiation into crust and "present mantle," while our results correspond, of 
course, only to the present mantle. Therefore, in order to perform a comparison with the results obtained in the 
previous Section, for each model we remove from the PM the crustal Th and U masses (which, in our global crust 
model of Sec. Ill, amount to Mxh = 15.30 ± 1.77 and Mu = 3.45 ± 0.57 at la, in units of lO^^ kg). 

The distribution of the remainder Th and U masses within the present mantle volume is subject to a lively de- 
bate, opposing homogeneous versus inhomogeneous (e.g., layered) models jH, For the purposes of this work, we 
consider two representative extreme scenarios, yielding "high" and "low" mantle geoneutrino rates. The "high rate" 
(homogeneous) scenario is obtained by subtraction of the Th and U crustal masses at the lower end of their la range, 
and distributing the remainder in the whole mantle at constant density. The "low rate" (inhomogeneous) scenario is 
obtained by subtracting from the PM the Th and U crustal masses at the upper end of their la range, and placing all 
the remainder in the so-called D" layer (250 km thickness) just above the core-mantle boundary. In both cases, aver- 
aged oscillations are included. Note that the various models are based on different assumptions or input values about 
the primitive chondritic material, which lead to further differences in the Th and U contents and in the associated 
radiogenic heat in the present mantle. 

Table V summarizes our estimated "low" and "high" Th and U mantle geoneutrino rates as derived from each mantle 
model, together with the associated total heat H{Th + U). Note that two models (Allegre et al. 1995, McDonough 
& Sun 1995) are practically identical for our purposes. In one model (Javoy et al. 2010), characterized by rather low 
radiogenic abundances, the Th content of the PM is slightly lower than the corresponding content of the crust in our 
"low" scenario; crustal subtraction would then yield slightly negative results, which have been just set to zero. Since 
not all models in Table V are endowed with PM error estimates, we do not quote individual errors for the present 
mantle rates; their spread is, however, indicative of the large theoretical uncertainties typically associated to mantle 
observables. 



TABLE V: Geoneutrino event rates derived from various models of the primitive mantle (PM), under different assumptions about tire Tli and U 
distributions in tire present mantle, leading to "low" and "high" rates. The first three columns characterize the original PM model in terms of Th 
and U masses. After crustal subtraction and redistribution of the remaining Th and U masses in the present mantle, we derive the oscillated Th 
and U mantle event rates, the Th+U heat and the Th/U ratio as reported in the last eight columns, for the "low" and "high" scenarios. 



Primitive mantle characteristics 




Present mantle 


, "low" scenario 


Present mantle. 


"high" scenario 


Model 


A/Th 


A/u 


i?.(Th) 


i?.(U) 


H(Th + U) 


Th/U 


i?.(Th) 


i?,(U) 


H(Th + U) 


Th/U 




[lO^^kg] 


[10^'kg] 


[TNU] 


[TNU] 


[TW] 




[TNU] 


[TNU] 


[TW] 




Turcotte & Schubert 2002 


3.62 


0.90 


2.7 


9.8 


17.0 


3.9 


3.9 


14.7 


19.0 


3.8 


Anderson 2007 


3.13 


0.78 


2.3 


8.4 


14.5 
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FIG. 4: Comparison of experimental constraints and model predictions in the plane charted by the Th and U mantle rates. 
Each model leads to extremal case of "low" and "high" rates, connected by lines to guide the eye. The KL+BX constraints 
are shown as na contours in steps of 0.5a. See the text for details. 

Figure 4 shows the same KL+BX mantle rate constraints as in the lower panel of Fig. 3, but with superimposed 
model predictions from Table V. In order to guide the eye, the "low" and "high" predictions for each model are 
connected by a straight line; moreover, the KL+BX na isolines are shown in steps of 0.5(7. It appears that the data 
prefer mantle models with relatively high radiogenic contents (e.g., Turcotte & Schubert 2002) and disfavor at ~ 2a 
those with low contents (e.g., Javoy et al. 2010). This indication is rather interesting and deserves to be checked and 
investigated with further data, given its potential impact on mantle model building. 

The near alignment of all model predictions in Fig. 4 reflects the narrowness of the Th/U mantle ratio in Table V, 
Th/U 6 [1.7, 3.9]. In this context, it makes sense to marginalize away the Th/U ratio within such range, and to 
express the resulting rates in terms of the total (Th+U) event rate. We find at la that 

i?(Th + U, mantle) ~ 23 ± 10 TNU (KL + BX, for Th/U £ [1.7, 3.9]) , (7) 

with nearly gaussian error distribution (not shown). At present, this represents our best estimate for the mantle 
geoneutrino flux, as derived by using inputs from particle physics (KL, BX, and oscillation data) and from Earth 
sciences (crustal data and mantle Th/U ratio). 

Finally, Fig. 5 shows a comparison between theory and data in terms of the Th+U mantle rate (in TNU) and 
radiogenic heat (in TW). The various model predictions, shown as lines connecting the "low" and "high" cases in 
Table V, can be compared to the mantle rate inferred in Eq. (O and shown as a horizontal la band. Notice that the 
la experimental error is already comparable to the spread of the theoretical rate predictions: therefore, improvements 
in the experimental accuracy by a factor of two or more, would allow a statistically significant model discrimination. 
At present we note that, at la, the data favor models with relatively high mantle heat iJ(Th + U), such as Anderson 
2007 and Turcotte & Schubert 2002. However, no model can be really excluded at ^ 2a. 

Signal- heat correlation plots as in Fig. 5 have been suggested earlier [l| as a way to constrain the radiogenic heat, 
whose preferred range should be determined by the intersection of the slanted "theoretical" band (i.e., the envelope of 
all models) with the horizontal "experimental" band. At present, given the large uncertainties affecting the bands in 
Fig. 5, we refrain from deriving a preferred range for the mantle heat — a task which is left to future studies involving 
more accurate data and, possibly, more extensive model surveys. We just observe that Fig. 5 suggests a lower bound 
for the Th+U mantle heat, which appears to exceed ~ 13 TW at la; then, adding no less than ~ 6 TW from Th+U 
crustal heat (see Sec. Ill A), we infer a tentative lower bound of ^ 19 TW (at la) for the total Th+U radiogenic heat 
in the Earth, compatible with the estimate '-^20 + 9 TW of [2] . 
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FIG. 5: Comparison of KL+BX constraints {la horizontal band) and model predictions (slanted lines) in the plane charted by 
the Th+U geoneutrino rate and radiogenic heat for the mantle. 



VI. SUMMARY AND CONCLUSIONS 



In this work, we have analyzed in detail the experimental total rates of Th and U geoneutrino events in KamLAND 
(KL) and Borexino (BX), including neutrino oscillation effects. We have calculated the crustal flux at the two detector 
sites, using updated information about the global and local Th and U distribution. After subtraction of the estimated 
crustal component from the total fluxes, we find hints for residual mantle components at > 1.5cr in both KL and 
BX. In the KL+BX combination, the statistical significance of the mantle signal reaches the 2.3a level. In particular, 
for typical Th/U mantle ratios, we estimate a total mantle rate of i?(Th + U) ~ 23 ± 10 TNU (including oscillation 
effects). The ±10 TNU error is comparable to the spread of rate predictions derived from various published models 
of the mantle. Among these, a preference is found for models with relatively high radiogenic contents (corresponding 
to present mantle Th+U heat > 13 TW at la). However, no model can be excluded at ^ 2a level yet. 

The accuracy of the results can be improved in part by further KL and BX data, and especially by prospective data 
from future experiments at different sites. If Th and U rates from several different detectors were available, one could 
estimate and subtract the crust components in all of them, in order to infer the corresponding mantle components. 
Should all mantle rates be the same within errors, the hypothesis of a common, isotropic mantle geoneutrino flux 
would be corroborated, and its value could be obtained by combining the results from all the experiments. Conversely, 
one should consider the possibility of an anisotropic mantle flux, or of incorrect estimates of (some of) the subtracted 
crustal fluxes. Mantle-dominated measurements at oceanic sites would provide crucial tests of the various options 
and additional constraints for mantle models. In any case, we shall learn a lot more about the Earth mantle from 
future geoneutrino data, supplemented by detailed descriptions of the crust at global and local scales, in a truly 
interdisciplinary approach. 
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Appendix A: Central values and uncertainties of estimated crustal rates 

In this work, a careful evaluation of crustal geoneutrino rates has been carried out (Sec. Ill) in order to infer the 
mantle signal by subtraction (Sec. IV). In this Appendix, we discuss possible effects or additional constraints that 
might slightly change either the central values or the uncertainties of our estimated crustal fluxes; we argue that such 
changes are not expected to significantly weaken the inferred mantle signal. 

1. Central values 

In Sec. Ill, crustal rates have been estimated by summing the LOG (local) component, based on a detailed description 
of the crust near the detector site, and the ROC (rest-of-the-crust) component, based on state-of-the-art information 
about the global structure and radiogenic contents of the crust. 

Of course, the global model of the crust used in Sec. Ill can be improved by adding new data or constraints, which 
might also lead to variations in the central ROC values. In particular, heat-flow measurements at local, regional, and 
global scales appear to provide particularly promising and independent constraints [Bj. In this context, we note that 
a recent heat-flow evaluation of the total radiogenic heat of the crust suggests a value _ff (Th -I- U-f K) ~ 7.5 ± 0.7 TW 
[S^ l (as quoted in [Sj), which is ^ 20% lower than our estimate 8.9 ±1.2 TW in Sec. Ill, although compatible within 
the quoted la errors. At present, it is not obvious if and how the above heat-flow estimate should be added as an 
additional input; in any case, its possible inclusion would presumably lead to an 0(10%) decrease of the adopted 
global Th and U crust abundances, and to an associated decrease of the estimated ROC fluxes in both KL and BX. 
The mantle flux estimated by subtraction would then slightly increase, and the hints for a mantle signal discussed in 
Sec IV would be strenghtened. 

Another issue concerns the approximation of averaged oscillations. While this approximation is certainly justified 
for the ROC component, it may be slightly inaccurate for the LOC component. We have re-calculated the LOC rates 
for the unaveraged oscillation probability in Eq. ([3]), and we find the following fractional variations with respect to 
the LOC rates in Table IV: +1% (Th) and -^3% (U) for KL, and -7% and -6% for BX. [In BX, a larger fraction 
of local fiux is suppressed by the first dip of the oscillation probability.] In terms of total (ROC-I-LOC) rates, all 
such variations are smaller than ±2.5% and can be ignored in practice (see also below). Note that, if the effects of 
unaveraged oscillations were statistically relevant for LOC estimates, then the analysis of geoneutrino energy spectra 
should also account for small L/E effects (not included in the present work). 

Based on the above considerations and on our educated guess, we surmise that future improvements in the local 
and global description of the crust are unlikely to change our estimated crustal rates by much more than 0(10%), 
which is also the size of the uncertainties in Table IV. Uncertainties of this (and even twice as large) size are negligible 
in the context of the present work, as we demonstrate below. 

2. Uncertainties 

In Sec IV we have inferred the mantle geoneutrino rates by subtracting the crust component from the total exper- 
imental rates, whose errors have been claimed to dominate the results. The irrelevance of the theoretical crust rate 
errors can be easily demonstrated in the case of KL, where the na experimental contours in the upper panel of Fig. 2 
are nearly elliptical and equally spaced (i.e., the correlated errors can be approximated by a bivariate gaussian). 

In order to simplify the argument and the notation, let us call (x, y) the total KL experimental U and Th rates 
shown in the upper panel of Fig. 2, in TNU units. The allowed contours are well approximated by ellipses centered 
at (27.8, 14.5) and with la error matrix 

\y pa^ay a^ 

where ax = 18.6, ay ~ 11.4, and p = —0.75. If we simply subtract the crustal rates (xc, j/c) — (20.2, 6.3) in 
Table IV with no errors, we get an estimated mantle signal {xm, Um) = (7.6, 8.2), with the same error matrix as 
above. The corresponding function [l^ equals 2.57 ~ 1.6^ at the point (0, 0), namely, the null hypothesis of "no 
mantle signal" is rejected at 1.6a in this gaussian approximation for the KL experimental errors (consistently with 
the non-approximated value of 1.7cr reported in Sec. IV). 

If the KL crustal rate errors in Table III are included with a (presumably positive) correlation pc in a crust error 
matrix (Tc, the mantle error matrix is augmented as cr^ = cr^ + cr^. The null hypothesis is then rejected at 1.58- 
1.59(7, depending on p G [0, 1]. Therefore, the inclusion of crustal rate errors changes the statistical significance of 
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the mantle signal in KL by a negligible fraction of one standard deviation (0.02(7 or less). If we repeat the above 
exercise with crustal errors doubled (or tripled) for the sake of conservativeness, then the significance of the mantle 
signal in KL is lowered by less than O.OSfi (or O.lSfi). For BX (middle panel of Fig. 3), the effect of crustal errors can 
only be smaller, since the experimental errors are significantly larger than in KL. Therefore, we surmise that the 2.4(t 
statistical significance of the mantle signal (Sec. IV) could be lowered by about 0.1-0.2cr in the worst cases. 

Smumarizing, variations or uncertainties of the estimated crustal rates at the level of 0(10%) can be practically 
ignored in the context of this work, whose results are dominated by the — much larger — experimental errors on geoneu- 
trino event rates. In the future, when theoretical and experimental errors will be comparable, a proper convolution of 
their distributions will be required in order to obtain reliable error estimates for the inferred mantle geoneutrino rate. 

Appendix B: Side results of the KL and BX data analysis 

In this Appendix we supplement the data analysis reported in Sec. II and in Fig. 1 with additional results, concerning 
total rates only, with no separation between mantle and crust components. As in [l3|, we re-express the experimental 
constraints in terms of total Th+U rates and Th/U abundance ratios, instead of separate Th and U geoneutrino rates. 
This alternative formulation is particularly useful to evaluate the effects of assuming, in first approximation, the same 
average Th/U ratio in both KL and BX [12|. The results are shown in the various panels of Fig. 5. 

Figure 6(a) shows the best fits and la contours in the plane charted by the total rate i?(Th-|-U) and by the average 
Th/U abundance ratio. In the upper panel, these observables are left free for both KL and BX. Both experiments 
provide similar (although weak) upper bounds on Th/U, but only KL sets a lower limit at present, Th/U ^ 1. In the 
lower panel, under the assumption of a common Th/U ratio in KL and BX, we obtain a best fit Th/U ~ 5, close to 
the typical chondritic value Th/U ~ 3.9 11], although with very large uncertainties (a factor of about four). 

Figure 6(b) shows the constraints on the Th/U abundance ratio in terms of standard deviations (the total rate 
being marginalized away). As in Fig. 6(a), the upper panel refers to the unconstrained fit, while the lower panel refers 
to the case with the same Th/U in KL and BX. It can be seen that, in combination, KL+BX can set for the first 
time a 2a upper bound on Th/U, although it is still very weak (Th/U ^ 10^). In this context, significant progress 
will require geoneutrino energy spectra with high statistics and low systematics, in order to better discriminate the 
different Th and U components and to constrain their ratio with higher accuracy. 

Finally, Fig. 6(c) shows the constraints on the total rate i?(Th + U) in terms of standard deviations No- (the Th/U 
ratio being marginalized away). The upper panel coincides with Fig. 1. The lower panel is rather similar, implying 
that current constraints on the total Th+U rate are not very sensitive to specific assumptions on the Th/U ratio. 




10 10" 10" 1 10 10' 10" " ™ °" 0" 
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(a) Best fits and Icr contours in the plane of (b) Constraints on the Th/U ratio in terms (c) Constraints on the Th+U rate in terms 
the Th+U rate versus the Th/U ratio. of standard deviations Na from the best fit. of standard deviations N„- from the best fit. 

FIG. 6: Analysis of KL and BX data in terms of the total Th+U event rate and of the average Th/U abundance ratio. In all 
the subfigures, the lower (upper) panel does (not) include the assumption of equal Th/U in KL and BX. 
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